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FOREWORD 


The  literature  indicates  that  the  use  of  microvave  energy  can  reduce 
freeze-drying  cycles  to  l/2  -  l/lOth  of  the  time  required  for  conventional 
(conductive,  radiant)  freeze-drying.  Two  of  the  major  drawbacks  to  design 
and  application  of  microwave  energy  to  freeze-drying  are  determination  of 
the  electric  field  available  for  dielectric  heating  and  corona  breakdown 
(gas  plasma  formation).  Corona  breakdown  can  cause  undesirable  effects  in 
food  products  during  freeze-dehydration,  such  as  deterioration  of  flavor 
components,  degrading  of  structure  and,  in  the  extreme,  burning  of  the  food. 

The  author  briefly  reviews  the  theory  of  microwave  gas  breakdown  and 
compares  theoretical  and  literature  breakdown  curves  for  air  and  noble 
gases  to  show  the  effects  of  pressure,  temperature,  frequency,  gas  composi¬ 
tion,  size  and  shape  of  cavity,  and  electric  field  strength.  The  author 
then  derives  an  equation  which  shows  the  effect  of  a  dielectric  load  in  the 
cavity  on  gas  breakdown,  and  shows  how  this  relates  to  heating  of  the  dielec¬ 
tric  load  and  to  electric  field  strength. 

Experimental  breakdown  curves  for  air,  water  and  carbon  dioxide  are 
given  and  related  to  theory  at  2450  MHz,  the  most  commonly  used  frequency 
for  food  materials.  (The  literature  has  very  sparse  data  on  air  breakdown 
at  2450  MHz  and  almost  none  on  water  or  carbon  dioxide  breakdown  at  ary 
frequency ) . 

Pressures  studied  were  in  the  range  of  0.1  to  20  Torr,  which  covers 
the  region  of  interest  in  practical  freeze -drying.  Temperatures  were 
essentially  ambient  (24°C).  A  single  cavity  approximately  two  wavelengths 
on  each  side  was  used  for  three  gases.  Amperex  DX206  magnetron  was  used 
with  a  calibrated  H- tuner  and  dummy  load  to  attentuate  the  power.  A  Pyrex 
vacuum  flask  inside  the  cavity  contained  the  corona,  and  gas  bleed  in  and 
hence  pressure  was  controlled  by  needle  valves  in  series.  A  bi-driectional 
coupler  sampled  forward  and  reflected  power  for  scope  display  and  calcula¬ 
tion  of  absorbed  power. 

The  agreement  between  theory  and  experiment  is  good,  both  for  air 
breakdown  and  the  effect  of  a  dielectric  load. 

The  results  obtained  show  for  an  arbitrary  microwave  freeze-dryer 
with  an  arbitrary  operating  pressure  and  food  load,  just  how  much  power 
can  be  applied  and  absorbed  by  the  food  without  corona  breakdown*  This 
plus  a  knowledge  of  maximum  mass  transfer  rates  of  water  vapor  across  the 
dried  food  layer  is  expected  to  enable  theoretical  optimization  and  predic¬ 
tion  of  microwave  freeze -drying  rates. 

Much  of  this  work  will  be  reported  at  the  Fifth  Annual  Meeting  of  the 
International  Microwave  Power  Institute,  at  The  Hague,  Netherlands,  on 
October  7-9,  1970. 
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ABSTRACT 


Pood  Freeze-drying  cycles  of  1/2  to  1/10  of  the  time  required 
ty  conventional  methods  (radiant,  conductive)  appear  feasible  using 
microwave  power.  A  major  problem  is  corona  breakdown,  which  occurs 
most  readily  under  the  pressure  conditions  of  freeze-drying  and  causes 
deleterious  effects  on  the  food. 


The  effects  of  pressure,  temperature,  frequency,  gas  composition, 
size  and  shape  of  cavity,  electric  field  strength,  and  dielectric  load 
are  explained  and  Interrelated  in  the  region  of  Interest  of  freeze¬ 
drying,  using  breakdown  curves  from  the  literature  derived  by  the 
authors.  Experimental  breakdown  curves  for  air,  water  vapor, 
carbon  dioxide  at  2450  MHz  sure  presented  compared  with  theory. 
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INTRODUCTION 


With  the  Amy’s  emphasis  on  freeze  dried  rations  and  the  large  number  of 
units  involved,  any  reduction  in  unit  cost  can  give  significant  savings.  This 
study  vas  undertaken  as  a  preliminary  step  to  applying  one  potential  method  of 
reducing  cost  -  that  of  combining  microwaves  and  freeze  drying.  The  usefulness 
of  microwaves  in  freeze  drying  has  been  suggested  and  reported  by  numerous 
people  1,2,3,^.  During  the  middle  and  end  of  the  freeze  dry  cycle,  microwaves 
by-pass  the  problem  of  heat  conduction  across  the  dried  food  layer  by  giving 
essentially  volumetric  heating  of  the  receding  ice  layer.  Thus  the  limiting 
step  in  the  process  becomes  mass  transfer  rathe  r  than  heat  transfer,  and  de¬ 
pending  upon  the  diffusivity  of  the  food,  the  cycle  time  can  be  reduced  by  a 
factor  of  two  to  ten. 

The  basic  drawbacks  of  the  process  have  also  been  reported.  Aside  from 
equipment  and  electricity  cost,  they  are: 

1.  Uneven  electric  fields  -  the  electric  field  is  not  uniform  for  any 
general  cavity.  The  use  of  mode  stirrers,  crossed  modes,  multiple 
inputs,  and  continuous  processing  can  reduce  this  problem  and  promote 
even  heating  of  the  frozen  food. 

2.  Meltback  -  if  microwave  energy  is  put  into  the  food  faster  than  the 
sublimation  -  mass  transfer  process  can  remove  it,  the  pressure  at 
the  ice  interface  will  rise  above  the  triple  point  and  the  ice  can 
melt.  Aside  from  interrupting  the  freeze  dry  process,  the  microwaves 
will  couple  selectively  into  the  water  rather  than  the  ice,  causing 


intense  local  heating,  accelerated  melting,  and  a  "runaway"  condition. 
Thus  the  dielectric  heating  rate  must  be  matched  to  the  rna-Hrmiw  allow¬ 
able  mass  flow  rate  to  optimize  the  process. 

3.  Corona  breakdown  -  at  reduced  pressures,  the  ionization  of  gas  molecules 
by  electrons  accelerated  in  the  electric  field  can  produce  a  plagma 
formation.  This  gas  plasma,  also  called  a  corona  or  glowball,  will 
consume  power  meant  for  food  heating,  affect  the  flavor,  and  can  even 
scorch  the  food  surface  if  sufficiently  intense.  The  breakdown  field 
strength  is  a  second  upper  limit  to  power  absorption  by  the  food 
(power  absorption  is  proportional  to  the  electric  field  squared). 

Both  2,  and  3,  are  related  to  pressure,  temperature,  electric  field 
intensity,  and  the  food  loads.  While  pressure,  temperature,  and  weight 
measurements  present  no  real  problem,  the  electric  field  is  a  difficult  quantity 
to  measure  directly .  The  applicable  theory  in  the  literature  does  not  relate 
well  to  multimode  cavities  and  high  power  absorption,  and  requires  rather 
sophisticated  measurements  of  VSWR  and  sweep  frequency  half  power  points,  methods 
not  very  applicable  for  industrial  monitoring  because  of  cost  and  untrained 
personnel. 

This  paper  will  instead  relate  food  loading  and  electric  field  to  power 
absorption,  which  can  be  measured  with  a  forward-reverse  power  meter  or  bi¬ 
directional  coupler.  From  this,  a  means  of  separating  food  power  absorption 
from  skin  loss  of  the  cavity  and  dielectric  loss  to  racks,  supports,  etc.  in 
the  cavity  will  be  shown.  Finally,  the  electric  breakdown  fields  will  be 
compared  with  theory  and  experiment  for  empty  cavities  and  loaded  cavities,  for 
varying  frequencies,  dielectric  constants,  pressures,  cavities,  and  gases. 
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Thus  we  tell  have  all  the  data  needed  for  process  optimization  except 
the  maximum  mass  flow  rate.  The  integration  of  this  paper  with  "the  mass 
transfer  literature  is  left  for  later  work. 

BASIC  THEORY  OF  GAS  BREAKDOWN 

A  corona  discharge ,  or  gas  plasma,  is  a  gaseous  mixture  of  neutral  gas 
molecules,  ions,  and  electrons.  Its  most  obvious  characteristic  is  a  visible 
glow  within  the  microwave  cavity  of  a  color  characteristic  of  the  gas  being 
used.  Air  gives  a  yellow-red  discharge,  water  a  reddish  glow,  and  carbon 
dioxide  gives  a  blue  color.  Figure  la  shows  a  relatively  dense  plasma  corre¬ 
sponding  to  a  high  power  level.  As  the  power  is  decreased  at  constant  pressure, 
the  plasma  weakens  as  in  lb  and  finally  reveals  the  "hot  spots"  of  maximum 
electric  field  as  in  Figure  lc  before  extinguishing. 

The  basic  concept  of  plasma  formation  or  gas  breakdown  is  simple:  free 
electrons  are  accelerated  by  the  oscillating  electric  field  according  to  the 
equation. 

a=accel.  cm/sec2 

1.  a=Eoe  sin  tut  Eo=max.  electric  field  volts/cm. 

m  e=electron  charge  coulombs 

m=electron  rest  mass  grams 
cu=field  frequency  radians/sec 

These  moving  electrons  can  have  elastic  or  inelastic  collisions  with  the 
gas  molecules.  In  elastic  collisions  the  electrons  bounce  off  the  atom  and 

gain  kinetic  energy,  but  no  change  occurs  in  the  internal  state  of  the  atom. 
After  numerous  elastic  collisions,  the  electron  may  gain  enough  energy  for  an 
inelastic  exciting  collision.  An  exciting  collision  increases  the  internal 
energy  of  the  atom  at  the  expense  of  the  electron's  kinetic  energy,  and  the 
electron  bounces  off  at  a  slower  speed.  The  energy  is  reemitted  as  a  quantum 
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hv  when  the  excited  atom  relaxes  to  its  ground  state. 

If  the  electron's  kinetic  energy  equals  the  ionization  potential  of  the 
gas,  an  inelastic  ionizing  collision  may  occur  which  releases  two  electrons, 
the  original  one  plus  one  stripped  from  the  neutral  molecule  ,  Actually,  this 
ionization  does  not  occur  as  a  single  step.  The  exact  mechanisms  are  not  cer¬ 
tain,  but  the  complexity  can  be  seen  in  the  twelve  step  process  postulated  by 
Campbell5  for  air. 

(A  fourth  type  of  collision  is  called  superelastic,  and  occurs  when  the 
electron  acquires  some  of  the  energy  of  an  excited  atom.  This  de-excitation 
is  often  important  in  gas  breakdown,  especially  with  pulsed  power,  but  will  not 
be  considered  in  this  paper.) 

When  the  production  of  electrons  in  ionizing  collisions  exceeds  electron 
loss  by  diffusion,  recombination,  attachment  and  wall  loss,  a  condition  similar 
to  a  chain  reaction  exists,  and  breakdown  occurs.  Most  breakdown  theory  uses 
the  criterion  that  when  the  rate  of  electron  production  equals  the  rate  of 
electron  loss,  breakdown  initiates.  The  difficulty  in  the  theory  is  finding 
the  appropriate  rates  for  the  various  mechanisms. 

A  full  explanation  of  these  mechanisms  is  far  beyond  the  scope  of  this 
paper.  Numerous  books  5,6,7, 8, 9  deal  ^th  the  subject  in  detail,  one  of  the 
best  being  MacDonald's  ^  This  paper  will  give  only  a  qualitative  summary  of 
the  processes  and  indicate  the  regions  where  a  particular  electron  loss 
mechanism  is  dominant. 

COLLISIONS  AND  ELECTRON  PRODUCTION 

Integrating  the  electron  acceleration  equation  1  gives  the  velocity 

2.  Velocity  =^e  Eo  cos  cut 
m  u) 
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Note  that  the  velocity  expresses  the  electron  flow,  or  current  I. 

Since  the  voltage  equals  Eo  sin  cut,  the  current  and  voltage  are  90°  out  of 
phase  and  the  average  power  gained  by  the  electrons  equals  zero: 

3.  Power  =  VI  cos  9  9  =  phase  angle  between 

current  and  voltage 
cos  90°  =  0 

V  =  voltage 
I  «  current 

The  collisions  of  electrons  with  gas  molecules  give  a  finite  value 
to  cos  9  and  allow  energy  transfer  from  the  field  to  the  electrons.  This 
energy  transfer  is  logically  dependent  upon  the  collisional  frequency  \fc , 
which  is  the  number  of  collisions  one  electron  has  in  one  second.  The 
magnitude  of  vc  depends  upon  the  cross  sectional  area  of  the  gas  molecule, 
the  electron  energy,  and  the  particle  density  of  the  gas. 

Temperature  and  pressure  affect  the  particle  density  and  hence  the 
collisional  frequency  as  well,  as  shown  by  rearranging  the  gas  law: 

4.  PV  =  N  ET  -  nET  N  =  number  of  molecules 

No  No=  Avogadro's  number 

6 . 02X10  ^  mole cule s/mole 

-5.  N  =  PNo  =  Particle  density 
V  ET 

Thus  increasing  pressure  or  decreasing  temperature  increases  collisional 
frequency.  Normally  all  pressures  for  gas  breakdown  are  converted  to  a 
standard  temperature,  either  0°C  or  20 °C. 

Assuming  vc  is  independent  of  electron  energy  (an  assumption  that  must 
be  used  with  considerable  caution)  an  average  value  for  air  can  be  expressed 
as : 
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6.  vc  =  5.3  X  lO9^  p  p  in  mm.  Hg  @  20°C. 

Obviously  the  collision  rate  will  interrupt  the  ordered  motion  of  the 
electrons  as  expressed  in  equation  2.  If  Vc  »u),  there  are  many  collisions 
per  field  oscillation  and  electron  velocities  can  never  reach  the  maximum. 

If  Vc«u),  there  are  many  field  oscillations  per  collision  and  the  velocities 
will  approach  those  expressed  by  equation  2. 

Correcting  for  collisions,  the  average  electron  drift  velocity  is 
expressed  by  : 


The  quantity  Ee  is  called  the  effective  field  and  can  be  considered  the 
equivalent  d.c.  field  for  fairly  high  collision  rates.  Experimental  drift 
velocities  have  also' been  reported  13>14>15. 

The  relation  of  vc  and  cju  gives  a  rough  idea  of  the  shape  and  location  of 
the  breakdown  curves.  Breakdown  is  easiest  when  the  frequency  of  the 
applied  electric  field  equals  the  collisional  frequency,  that  is  when  vc  =  tu. 
When  vc  (and  pressure  as  expressed  by  equation  6)  is  above  or  below  that  point, 
the  average  electron  energy  on  collision  is.  reduced  and  a  higher  electric  field 
is  needed  for  breakdown.  We  expect  a  plot  of  breakdown  electric  fields  versus 
pressure  to  be  concave  upwards,  with  the  minimum  at  vc  =  <ju,  or  expressed 
differently  for  air : 

8.  p  =  2tt  f 

5.3  X  10s  f  in  Hertz 
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However,  instead  of  collisional  frequency  we  need  expressions  for  the 
ionization  rate  to  balance  against  electron  loss  to  satisfy  the  breakdown 
criterion  that  rate  of  electron  production  equals  rate  of  electron  loss. 
The  ionization  rate  can  be  related  to  either  collisional  frequency  vc  or 
drift  velocity  Vd: 


9.  vi  =  hi  vc 
10.  v,i  =tfv 


,  _  probability  of 

*  ionization 


Y  —  ion  pairs  produced 
cm  travel 


Values  for  the  constants  are  given  in  the  literature  -^A7  >18,19, 20  anc^ 
show  an  energy  dependence  which  further  complicates  detailed  analysis.  To 
get  a  feel' for  the  numbers  involved,  an  hi  for  hydrogen  with  electron  energy 
of  18  electron  volts  is  given  as  .01.  vc  for  hydrogen  is  approximately 
5.9  X  lO^Cp,  so  at  1  mm  Hg  Vj_  =  5.3  X  109X  10a=  5.3  X  107  ionizations/sec. 


ELECTRON  LOSS  MECHANISMS 

Recombination 

If  an  electron  collides  with  a  positive  ion,  they  may  recombine  to  form 
a  neutral  molecule.  Generally  the  positive  ion  concentration  before  breakdown 
is  too  low  for  the  recombination  rate  to  be  significant.  It  may  be  important, 
however,  in  determining  decay  times  for  pulsed  breakdown.  The  rate  vr  is  given 
by: 


11.  vr  -  hp  vc 


hr  ~  probability  of  recombination 

per  collision 
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ATTACHMENT 

An  electron  can  also  attach  to  a  neutral  atom  or  molecule  and  be  lost  to 
the  breakdown  process.  Even  though  the  negative  charge  of  the  ion  is  equal  to 
that  of  an  electron,  the  mass  is  so  much  greater  that  the  drift  velocity  given 
by  equation  7  is  neglibible. 

The  attachment  rate,  like  the  ionization  rate,  has  been  related  to  both 
collision  frequency  and  drift  velocity: 

12.  v*  =  hftvc  =  efficiency  of  attachment 

13.  v.  =  ccy  oc  =  number  of  attachments  per 

^  cm  travel 

High  attachment  rates  will  be.  observed  with  molecules  whose  outer 
shells  are  nearly  filled,  such  as  chlorine  or  oxygen.  Halogenated  fLuorocarbons 
display  the  same  effect.  Noble  gases  with  their  stable  octet  have  very  small 
attachment  rates.  The  electron  affinity  energy  is  a  good  measure  of  relative 
attachment  rates  of  gases  for  which  oC  or  h_  are  unknown. 

Attachment  is  important  in  air  breakdown  because  of  oxygen  21, 22, 23^ 
is  the  dominant  loss  mechanism  only  at  ’’high"  pressures  and  correspondingly 
"high"  values  of  vc>  A  pressure  is  considered  "high"  when 

14.  p  A  >  10cm-mmHg  A  in  cm;  p  in  mm  Hg 

(A,  the  characteristic  diffusion  length,  is  a  measure 
of  the  cavity  size  and  shape  and  will  be  discussed  in 
the  section  on  diffusion) 

When  the  inequality  14  holds,  there  are  so  many  collisions  per  field  oscillation 
that  energy  modulation  is  negligible  and  Ee/p  ~  32  volts/cm  -  mm  Hg  regardless 
of  frequency  for  air.  This  determines  the  right  hand  side  of  the  breakdown  curve. 
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AMPLITUDE  LOSS 


When  the  amplitude  of  electron  oscillation  equals  the  size  of  the  cavity, 
electrons  will  be  lost  to  the  conductive  walls  and  will  not  be  available  for 
ionization,  causing  higher  breakdown  fields.  At  very  "low"  pressures,  the 
electrons  may  strike  the  walls  with  sufficient  energy  to  cause  re-emission  of 
another  electron  to  the  ionization  process,  thereby  lowering  the  breakdown  field. 
The  left  end  of  the  breakdown  curve  is  expected,  then,  to  go  through  a  maximum 
and  decrease  again  at  very  low  pressures. 

Theoretical  calculations  of  when  amplitude  loss  becomes  dominant  is 
difficult  10 ,  but  empirically  it  would  be  expected  at  an  inflection  point  in 
the  curve  or  when  the  first  derivative  of  the  curve  eauals  zero.  Because  of 
the  large  cavities  employed  in  commercial  freeze  drying,  amplitude  loss  usually 
will  not  be  dominant  in  the  normal  pressure  range  of  freeze  drying. 

DIFFUSION: 

Since  attachment  is  dominant  at  "high"  pressures  and  amplitude  loss  is 
dominant  at  "low"  pressures,  diffusion  is  said  to  dominate  the  "middle"  range 
of  pressure,  which  is  the  region  of  interest  for  freeze  drying  (around  1  mm  Hg). 
Diffusion  of  electrons  should  be  visualized  as  a  movement  superimposed  on  the 
electron  oscillations  caused  by  the  electric  field.  If  the  centers  of  oscil¬ 
lation  are  traced,  they  will  move  in  the  direction  of  decreasing  electron 
concentration  (the  walls). 

If  we  take  the  electron  concentration  as 

15.  V>  e  =  Ne  Ne  =  number  of  electrons 

Vol  Vol=  volume,  cm3. 
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then  the  driving  force  of  diffusion  is  the  gradient  of  electron  concentration 

16.  V"  e  =  2S£  +  +  |SS  : 

bx  by  bz 

The  flow  of,  electrons  by  diffusion  is  expressed  then  by  Fick's  Law: 


17.  j  =-VD,( 


Q  =  number  of  electrons/sec. 

A  =  cross  sectional  area  for  electron  flux 
D  =  diffusion  constant 


Note  that  equation  17  is  similar  to  the  heat  conduction  formula 

18.  =  —  k  q  =  Bttt/hour 

^  A  =  area  for  heat  flux 

k  =  constant  of  thermal 
conductivity 
dt/dx  =  temp,  gradient 


The  difference  is  that  heat  conduction  is  energy  transfer  and  diffusion  is 
mass  transfer. 

The  rate  of  change  of  electron  concentration  by  diffusion  then  becomes 

19.  =.y?(Dne) 

or  for  D  independent  of  spatial  variation 

2°.  =  DV8V\e 


Continuity  Equation  (Electron  balance) 

Using  the  criterion  for  breakdown  that  the  electron  production  rate  equals 
the  sum  of  the  loss  rates,  we  can  combine  equations  9,  12  and  20  to  give 

6  jlO 

vine  —  vane  — .  o. . 


21.  5ne  = 
dt 
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v  and  v  have  the  Vie  term  because  they  are  both  rates  for  a  single  electron. 

vi  a 

Rewriting  21  we  get 

22.  (vi  -  va)|fle  = 

Unfortunately,  the  values  of  D  for  electrons  moving  in  high  electric 
fields  are  not  generally  known.  Theoretical  calculations  of  D  depend  on 
electric  field,  pressure,  and  the  electron  velocity  distribution  function 
(the  average  drift  velocity  Vd  is  not  sufficient).  MacDonald10  has  treated 
this  problem,  but  less  detailed  analysis  makes  the  assumption  of  the  ioniza¬ 
tion  rate  being  independent  of  spatial  variation  (i.e.,  the  electric  field 
is  uniform  throughout  the  region  of  interest). 

With  this  assumption,  both  sides  of  22  may  be  set  equal  to  a  constant 
-  such  that 

23.  V  3  e  -  -  tes. 

The  value  of  tf/D  depends  upon  the  size  and  geometry  of  the  cavity.  Brown7 
and  MacDonald10  give  the  following  relations  for  different  geometries: 

For  right  cylindrical  cavity  with  conductive  walls: 

2L  ^/D  =(2.405  \s  +  /n_ \  Rs  radius  of  cylinder 

V  R  /  '  X »  height  of  cylinder 

For  rectangular  parallelepipeds  in  the  dominant  mode: 

25.  Vd  =  nfX®  +  Tiff*  +  n^Z3  X,  Y,  Z  dimensions  of  sides 

Because  ^D  has  units  of  cm-3  a  characteristic  diffusion  length ,  A.  . 
is  defined  such  that 
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26.  JT/d  =  l/A3 

For  the  case  of  infinite  parallel  plates ,  the  diffusion  length  is  simply 

27.  A  —  d/n  d  —  plate  separation,  cm 

Most  micro-wave  breakdown  work  has  been  done  with  A  a  fraction  of  a 
centimeter  or  a  couple  of  centimeters  at  most,  in  cavities  designed  for  a 
single  mode,  most  often  the  TE01 Q34  .  But  the  general  case  for  an  industrial 
freeze  drier  is  a  large  multimode  cavity  with  a  geometric  A  much  larger.  How¬ 
ever,  with  a  multimode  cavity  the  assumption  of  uniform  electric  field  implicit 
in  equations  24  >  25  and  26  no  longer  holds. 

The  true A is  much  smaller  than  the  geometry  indicates,  because  if  electrons 
move  from  a  high  electric  field  region  to  a  lower  one  they  are  as  effectively 
lost  to  breakdown  as  if  they  diffused  to  the  walls.  The  electric  field  density 
as  expressed  by  the  dominant  mode  of  the  cavity  is  also  important  in  that  it 
affects  the  assumption  of  vi  being  independent  of  spatial  variation® 6 , 

A  for  large  cavities  cannot  normally  exceed  one  wavelength  of  the  applied 
field10,  since  this  represents  the  maximum  periodic  spatial  variation  of  the 
electric  field.  The  large- cavity- A can  be  much  smaller  if  localized  high 
intensity  fields  exist.  Such  localized  fields  cause  a  sphere  or  ovoid  of 
breakdown,  the  so  called  "glowball",  and  the  approximate  diffusion  length  is 
given  by 

28.  A.  =  r/n  r  =  radius  of  "glowball" 

Sharp  points  or  discontinuities  in  the  microwave  freeze  dry  cavity  will  cause 
such  localized  fields  and  must  be  avoided,  both  to  minimize  breakdown  and  to 

ensure  even  dielectric  heating  of  the  food. 
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We  have  already  discussed  the  effects  of  temperature  and  pressure  and 
are  now  in  a  position  to  qualitatively  understand  the  effects  of  frequency, 
gas  composition,  and  cavity  size  and  shape  upon  breakdown  curves, 

GAS  COMPOSITION 

We  expect  different  gases  to  have  differing  breakdown  curves,  A  high 
ionization  potential  gas  requires  high  energies  and  correspondingly  high 
electric  fields  for  breakdown,  while  a  gas  like  neon  requires  relatively  low 
fields.  A  gas  with  a  high  electron  affinity  energy  such  as  chlorine  or  freon 
should  have  a  high  attachment  coefficient  and  also  require  high  fields,  while 
gases  like  argon  should  require  low  fields.  A  gas  with  a  large  collisional 
cross  sectional  area  should  have  a  relatively  high  collisional  frequency  vc  , 
and  the  minimum  breakdown  should  occur, at  relatively  lower  pressures,  as 
shown  by  equation  8. 

Figure  2  illustrates  these  points,  and  Table  1  gives  some  applicable 
physical  constants. 

Table  1  -  Selected  Constants  of  Gases 
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Effective  collision 

First  loniza- 

Gas 

Molecular  diameter-viscosity 

Molec.  Dia, 

tion  Potential 

o 

o 

°2 

2.96  A 

4.3  A 

12.1  e.v. 

0 

0 

N2 

3.16  A 

5.5  A 

15.6  e.v. 

0 

o 

H2 

2.18  A 

3.1  A 

15.44  e.v. 

o 

0 

Ar 

2.86  A 

3.7  A 

15.77  e.v. 

o 

o 

ro 

-F- 

8 
>  o 

0 

7.6  A 
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KLeeding  in  a  gas  with  higher  breakdown  strength  than  air  is  one 
method  of  increasing  the  allowable  electric  field  for  microwave  freeze 
drying.  Note  that  a  gas  may  have  a  breakdown  field  higher  than  a  second 
gas  in  one  pressure  region  and  lower  in  another  region.  Care  must  be  taken 
to  do  comparisons  in  the  pressure  range  of  interest. 

Unfortunately  breakdown  curves  are  not  known  for  very  many  gases.  A 
large  part  of  the  literature  concerns  itself  with  low  breakdown  strength 
noble  gases  because  of  the  good  agreement  between  theory  and  eiqieriment. 
Even  very  sophisticated  theory  fails  in  predicting  breakdown  curves  for 
polyatomic  molecules  because  of  computational  difficulties  in  establishing 
the  electron  velocity-energy  distribution  function,  and  because  of  the  lack 
of  accurate  physical  constants  (h  ,  h  ,  h  ,  v  ,  etc.).  Thus  a  realistic 
approach  to  choosing  a  bleed— in  gas  would  be  trial  and  error,  guided  by 
the  qualitative  principles  stated  above. 

FREQUENCY 

Since  the  minimum  in  the  breakdown  curve  occurs  approximately  where 
(l)  ~  Vgj  we  expect  the  curve  to  shift  to  the  right  for  higher  frequencies. 

A  second  result  of  higher  frequencies  is  an  upward  shift  of  the  curve, 
because  of  the  shorter  half  cycle  time  for  acceleration  of  the  electrons. 
Recalling  the  velocity  equation 

2.  Velocity  =  ^  ^  cos  <ut 

we  see  that  increasing  the  frequency  (l)  decreases  the  magnitude  of  velocity 
and  energy  modulation,  thereby  decreasing  the  rate  of  electron  production. 
Figure  3  shows  the  curve  shift  due  to  frequency. 
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Working  at  higher  frequencies  would  yield  a  double  advantage  of  higher 
breakdown  fields  and  greater  heating  of  the  dielectric  according  to  the  familiar 
power  equation. 

.14 

29.  P  =  8.85  x  10  ■  2  tt  f  e"  r  E? 

P  =  power  watts/  cm  3 
E  =  electric  field  volts/cm 
e"r  =  relative  loss  factor 
f  =  frequency  Hertz 

Higher  frequencies  do  have  the  drawback  of  smaller  half  depth  penetration; 
but  for  the  case  of  frozen  foods  as  in  freeze  drying,  the  penetration  at  2450  MHZ. 
is  still  deep  enough  to  assume  as  good  volumetric  heating  as  at  915*^  MHZ. 

The  loss  factor  e"  r  is  also  a  function  of  frequency,  but  for  most  foods 

the  factor  fe"r  increases  for  increasing  frequency  28 ,  so  the  power  absorption 

without 

is  still  increased  or  remains  constant/ changing  the  electric  field.  When  we 
also  square  the  allowable  breakdown  fields  shown  in  figure  11  for  air  at  915 
and  2450  MHZ  and  insert  into  equation  29,  the  advantages  of  higher  frequencies 
become  apparent. 

The  conclusions  of  various  researchers  working  at  915  MHZ.  that  very  low 
operating  pressures  are  needed  12 and  that  the  economics  of  microwave  freeze 
drying  are  unfavorable  28  do  not  necessarily  hold  at  2450  MHZ,  the  other  common 
ISM  frequency  for  food  processing.  For  example,  at  1  mm  Hg  and  for  cm) 
the  breakdown  field  for  2450  MHZ  is  about  230  volts/cm  and  for  915  about  100  volts/ 
cm.  The  ratio  of  maximum  power  absorption  without  breakdown  is  therefore: 

30.  P— 24.50  -/2i50  MHZ\  /e»r  2450\  (E  breakdown  2450)3 

P  915  915  MHZ  j  ye"r  915y  (E  breakdown  915 )2 
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Taking  representative  eMr  for  ground  beef  for  1,000  and  3>000  MHZ  89 , 

the  ratio  P  ?450  approximately  equals  5.3. 

P  915 

Clearly  frequency  has  a  strong  effect  on  economics,  and  2450  MHZ  may  be 
more  suitable  than  915  MHZ  for  freeze  drying. 

Effect  of  Cavity  Size 

For  a  given  frequency  and  gas,  we  expect  the  minimums  of  the  curves  for 
different  diffusion  lengths  (A)  to  fall  at  about  the  same  pressure.  The  electric 
fields  for  "high"  pressures  are  the  same,  since  for  P  A  >  10  attachment  is  the 
dominant  loss  mechanism  and  breakdown  fields  are  relatively  independent  of 
cavity  size.  At  "low"  pressures  where  the  amplitude  of  oscillation  can  reach 
its  maximum  (because  for  vc«cu  there  are  many  field  oscillations  per  collision) 
we  expect  lower  breakdown  fields  for  a  larger  A.  All  this  means  is  that  for  a 
larger  A  there  is  a  larger  region  of  uniform  intense  field  to  give  greater 
acceleration  and  energy  to  the  electrons,  and  that  the  cavity  is  larger  and  can 
contain  these  larger  oscillations  without  amplitude  loss  to  the  walls.  Figure 
4  shows  these  effects. 

A  increases  beyond  one  centimeter,  the  differences  in  the  curves 
become  less  pronounced.  This  is  helpful  because  for  large  cavities  A  will  be 
large  and  difficult  to  determine  precisely  from  estimation  of  the  radius  of  the 
"glowball".  An  alternate  method  is  to  estimate  the  point  where  the  breakdown 
curve  becomes  linear  and  assume  A  =  10/p.  Both  methods  are  crude  but  usable 
when  a  geometric  A  fails. 
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Summary  Simplified  Breakdown  Theory 

If  vc  was  independent  of  electron  energy,  all  breakdown  data  for  a  gas 
could  be  presented  on  a  two  dimensional  graph  using  the  lumped  parameters 
Ee  A  vs.  p  A,  as  in  Figure  5.  (Remember  that  Ee  includes  a  frequency  term, 
so  that  electric  field,  frequency,  pressure,  and  cavity  size  are  all  represented.) 
Figure  5  works  very  well  for  the  conditions  of  p  A  >  10,  when  Ee/p  —  32  and  for 
p  ^  >  100,  when  Ee  =  E  and  breakdown  becomes  independent  of  frequency.  At  lower 
values  of  p  A  the  assumptions  no  longer  hold  and  vc  becomes  energy  dependent. 

A  more  complex  semi-empirical  system  by  MacDonald  10relates  all  available  break¬ 
down  data  for  air.  Again,  the  theory  does  not  apply  to  breakdown  in  other  gases 
and  must  be  used  cautiously  in  predicting  freeze-dry  breakdown  fields  because 
of  the  water  vapor  partial  pressure. 

Lastly,  none  of  the  literature  theory  or  experiment  allows  for  a  di¬ 
electric  load  in  the  cavity.  The  author's  work  in  this  paper  will  therefore 
derive  and  test  expressions  for  the  effects  of  cavity  loading  on  breakdown 
fields.  Also,  no-load  curves  for  air,  water  vapor,  and  C02  for  2450  MHZ  as 
determined  by  the  author  in  a  multimode  cavity  will  be  presented  and  compared 
with  the  qualitative  principles  mentioned  before. 

Experimental  Work 


2450  MHz  was  used  for  all  of  the  present  work  because  of  the 
assignment  of  this  frequency  to  I.S.M.  use,  and  the  availability  and  common 
use  of  2450  MHZ  sources  for  food  processing.  Figure  6  shows  the  equipment  in 
block  form. 

The  source  A  is  an  Amperex  DX-206  1  K W  magnetron.  It  is  connected  via 
1  5/8-inch  coax  B  to  a  waveguide  circulator  C  which  shunts  the  reflected 
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power  to  a  dummy  load  D  preventing  reflection  to  the  source.  An  H- tuner  E 
varies  the  admittance  of  a  second  dummy  load  F  to  the  forward  power,  giving 
variable  attenuation.  A  bidirectional  coupler  G  samples  forward  and  reflected 
power,  and  a  switch  H  sends  one  or  the  other  RF  signal  to  a  rectifier  I  which 
gives  a  voltage  output  for  display  on  an  oscilloscope  J. 

A  waveguide  switch  K  directs  the  forward  power  to  either  the  cavity  0 
or  a  water  load  L.  A  thermopile  M  gives  a  voltage  proportional  to  the  temper¬ 
ature  rise  in  the  water  load,  which  is  read  on  a  millivoltmeter  N.  Figure  7 
shows  the  overall  equipment  set  up. 

An  antenna  beams  the  power  into  the  cavity.  The  vacuum  system  is  a  Pyrex 
flask  P  placed  within  the  cavity,  with  gas  bleed-in  from  the  sources  R  con¬ 
trolled  by  needle  valves  in  series  Q.  A  Stokes-McLeod  gauge  S  reads  pressures 
of  0-5mm  and  a  Wallace  and  Tiernan  diaphragm  gauge  T  reads  pressures  of  9-100  mm. 
A  vapor  trap  U  cooled  by  dry  ice  and  ethanol  prevented  water  contamination  of 
the  vacuum  pump  V  during  water  breakdown  runs.  Figure  8  shows  the  vacuum  flask- 
cavity  set  up. 

The  equipment  used  has  limitations.  First,  as  noted  by  Towne  29  a  magnetron 
is  a  rather  noisy  source  when  speaking  of  spectroscopy  or  breakdown,  and  most  of 
the  literature  work  was  done  with  ultrastable  monochromatic  sources.  A  noise 
source  can  have  lower  breakdown  fields  than  a  pure  source,  but  only  if Hufetmo. 

For  commercial  sources  of  2450  MHZ  +  50  MHZ,  this  effect  can  be  ignored. 

A  second  equipment  limitation  is  that  with  constant  bleed-in  of  gas,  break¬ 
down  occiirs  in  a  gas  flow  rather  than  under  the  static  conditions  of  most  lit¬ 
erature  work.  Freeze  drying  of  course,  involves  gas  flow,  and  the  author  wished 
to  investigate  the  effects.  Skinner  and  Brady  30 related  flow  effects  to  the 
critical  velocity  of  the  gas  while  duplicate  runs  by  the  author  with  flow 
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and  stagnant  gas  showed  no  significant  difference.  Care  must  be  taken  in 
static  breakdown  however,  so  that  all  the  products  of  the  plasma  are  removed 
and  fresh  gas  added  before  the  next  run.  Failure  to  do  this  can  give  break¬ 
down  fields  up  to  10$  too  low31. 

Another  source  of  error  was  temperature  and  composition  control.  Air  was 
no  problem  since  dessicated  ambient  air  was  used.  The  vacuum  boiled  distilled 
water  used  as  a  water  vapor,  source  had  dissolved  gases  that  had  to  be  distilled 
off,  and  small  air  leaks  in  the  system  contaminated  both  the  water  and  carbon 
dioxide  runs.  The  temperatures  of  water  and  C02  gas  were  below  ambient  at  the 
sources  and  were  warmed  by  conductive-convective  heating  through  long  lead-in 
tubes  to  the  cavity,  but  this  method  is  far  from  precise.  Thus  while  the  air 
breakdown  curves  are  quantitative,  the  other  two  should  be  considered  as  giving 
qualitative  breakdown  relative  to  air. 

A  final  equipment  limitation  is  that  the  magnetron  power  supply  used  in  the 
author's  work  is  a  full  wave  rectifier  type  which  gives  pulsed  rf  output  at  twice 
the  line  frequency  (i.e.,  2  x  60  HZ  =  120  pulses/sec).  The  breakdown  fields  for 
pulsed  rf  are  greater  than  for  CW,  a  nd  Copson32  has  suggested  pulsed  power  as  a  way 
to  increase  the  electric  field  and  dielectric  heating  in  microwave  freeze  drying. 

The  breakdown  fields  for  pulsed  rf  are  greater  than  those  for  CW  because 
there  are  fewer  collisions  in  which  the  electron  can  acquire  ionization  energy. 

The  energy  acquired  in  one  pulse  is  degraded  by  elastic  collisions  before  the 
next  pulse  occurs,  so  the  electron  chain  reaction  must  initiate  in  a  single  pulse. 

The  theoretical  approach  to  pulsed  breakdown  is  difficult  because  the  continuity 
equation  must  be  solved  for  transient  rather  than  steady  states. 

The  experimental  work  on  pulsed  breakdown  has  been  summarized  by  MacDonald1'0 
as  in  figure  9,  which  uses  the  ratio  of  (E  pulsed  breakdown  field/ECW  breakdown  field) 
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as  a  function  of  cycles  per  pulse,  allowing  comparison  of  different  frequencies 

and  pulse  lengths.  As  might  be  expected,  long  pulses  give  values  approximating 
CW. 


Recalling  that  absorbed  power  is  proportional  to  the  electric  field  squared 
and  correcting  for  the  finite  time  between  pulses  when  power  absorption  equals 
zero,  we  can  write 


31.  -  pulse  (Eprms)3  tp  feprms)3 

P  cv  ( Ecwrms  )8  Tp  *  (EcwrmsT^ 


(tp) 


(tp) 


Where  tp  -  pulse  time  per  cycle  with  non  zero  power. 

Tp  =  total  pulse  time  peak  to  peak  =  — - - - 

pulse  repetition  rate 

=  1/tp 


Therefore  for  pulsed  power  to  be  an  advantage,  the  quantity  (P  pulse/P  cw) 
of  equation  31  must  be  greater  than  unity.  For  the  author's  work  tp/Tp  =  0.743 
(determined  from  oscilloscope  display  of  forward  power  pulses.)  Taking  that 
portion  of  the  pulse  within  10#  of  the  peak  voltage,  there  are  3.76  x  10s  cycles/ 
giving  Ep/E  cv  ~  1.04,  and  a  ratio.  Pp/Pcw  of  0.804.  Much  shorter  pulse 
lengths  are  thus  needed  to  make  pulsing  advantageous.  For  example,  with  1  ^  sec 
pulses  at  2450  MHZ,  250,000  pulse/sec  would  be  necessary  just  for  Pp  to  equal 
Pew.  Normal  pulse  units  do  not  go  over  100,000  pulses/sec,  and  the  cost  of 
microwave  sources  with  higher  pulse  rates  and  the  large  duty  cycles  involved 

would  tend  to  offset  any  increase  in  absorbed  power.  Pulsing  does  not  seem  to 
be  the  answer. 


Determining  Electric  Field  Strength 

We  still  need  a  method  of  relating  power  to  the  electric  field  strength. 

Most  literature  work  measures  the  VSWR  on  resonance,  the  cavity  Q  and  the  incident 
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power  at  the  instant  of  breakdown.  This  method  though,  is  only  applicable  to 
empty  cavities  operating  in  a  single  dominant  mode.  The  author  derived  instead 
a  less  precise  but  more  general  form  involving  absorbed  power  allowing  for  loads. 
Since  forward  and  reflected  power  can  be  measured  by  a  bidirectional  coupler  or 
power  meter,  assuming  negligible  VSWR  between  the  detector  and  the  cavity  we  can 
write 


32.  Power  in  -  Power  out  =  Power  absorbed  in  cavity 

33.  P  absorbed  =  P  skin  loss  +  P  dielectric  loading  loss 

Since  the  skin  loss  is  analogous  to  a  resistance  loss  in  an  LC  circuit,  we  postu¬ 
late  a  form  for  the  first  loss  of 


34.  "  P  skin  loss  =  KEP  where  K  = - - - —  =  admittance 

resistance 

The  value  of  K  can  be  determined  for  an  unknown  cavity  by  measuring  P  in  and 
P  out  for  the  empty  cavity  for  a  known  breakdown  field.  For  an  arbitrary 
multimode  cavity,  A  will  not  be  known,  so  an  arbitrarily  high  pressure  is  used 
to  make  p  A  >  10,  giving  Ee/p  =  32  volts/cm  mm  Hg.  If  p^>  100,  the  breakdown 
also  becomes  independent  of  frequency  and  available  literature  curves  for  other 
frequencies  may  be  used  directly. 

Rearranging  equation  34. ,  K  is  simply 

35.  K  =  P  in  -  P  out 

IP 


For  a  cavity  not  yet  built,  a  reasonable  value  of  K  can  be  calculated  from 


1 


P°Qu 

I'M, 


or 


Po  =  Eo2  HI2. 

Qu 


KEP 
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Which  means  that 


37.  K  =  M 

Qu 


Where: 


Wo  resonant  freq, ,  radians/sec 

Po  =  power  absorbed  at  resonance,  watts 

Qu  =  2  tt  energy  stored  in  cavity 
energy  lost  in  cavity 

38.  T|  =  b  ec  V 

b  =  constant  of  electric  field  distribution 
e0=  free  space  permittivity  8.84XL014  farad/cm 
V  =  volume  of  cavity  cm  3 


The  theoretical  Qu  is  given  by  Harvey  33  for  various  dominant  modes. 

Assuming  a  T  E  202  mode  (a  reasonable  set  of  eigen  values  for  the  author's 
cavity),  we  get 

39.  ABC  [fr/A)8  +  (n/c)2  ]3/b 

0-/A)2  C  (A+2B)  +  (t  C  +  2B 

Where:  m,  n,l,  are  eigenvalues  of  T  E  m  n 

A,  B,  C  are  dimensions  of  cavity  A=B=C=24. 5cm 
6s  =  skin  depth. 

\  _ 

n  ~  ^ree  space  wavelength  =  12.25  cm 

Substituting  appropriate  values  in  37,  38  and  39  and  taking  b  =  (0.707)2 
for  the  average  electric  field  squared  for  half  sinusoidal  distribution,  we 
get  a  theoretical  K  of  0.37  x  10“3  versus  an  empirical  K  of  O.39  x  1CT3 
determined  by  equation  37. 
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This  is  a  very  close  correlation,  considering  the  exact  £s  of  the  cavity  was 
not  known  and  that  only  the  dominant  cavity  mode  was  considered.  The  theore¬ 
tical  method  should  give  reasonable  values  of  skin  loss  for  theoretical  design 
of  microwave  freeze  dry  cavities. 

Large  cavities  have  a  real  advantage  in  reducing  skin  loss,  since  Qu  is 
roughly  proportional  to  the  volume  divided  by  the  surface  area.  Taking  our 
cubic  cavity  of  side  X,  the  volume  equals  X?  and  the  area  equals  6  X? .  Thus 
Qu  increases  with  X  as 

40.  X3--  X 

ZIP  6 

and  the  large  Qu's  associated  with  larger  X's  give  smaller  K's  according  to 
equation  37 ,  and  correspondingly  smaller  skin  losses  for  a  given  electric 
field  by  equation  34. 

As  for  the  dielectric  loss,  we  refer  again  to  the  power  equation  29. 

41.  P  =  2  tt  f  e"  r  8.85xlCrr14  2  Volume  dielectric 

Note  that  we  want  total  watts  loss,  so  we  have  multiplied  3}  by  the  volume 
of  dielectric.  We  have  also  introduced  a  subscript  E^,  because  the  field 
inside  the  dielectric  is  less  than  the  free  space  electric  field  Ef  by 

42.  Ed  - 

We  must  keep  dearly  In  mind  that  the  breakdown  and  skin  loss  are  determined  by 
Ef  and  the  dielectric  heating  by  Ed,  all  expressed  as  rms  values.  Also,  41. 
must  be  corrected  for  pulsing.  Combining  31. ,  41.  and  42.  we  get 
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-3  p  =  2  TT  f  ey  8.85^L0~14(Efn)s  Vol.  Diel.  tp 

(®"r)3  +  (e  r)2  ^  Tp 

Where  Ef  refers  to  the  rms  of  the  pulse  alone,  since  the  term  iE  corrects 
p  Tp 

for  the  zero  power  portion  of  the  cycle.  Combining  35,  4l,  and  42  we  get  for 
CW  input  to  a  cavity  with  n  dielectric  loads 


44.  Efcw  = 


Fin-  Pout 


n 

Z 


K  +  L  2  TT  f  eVi  8.8$3cL0~l4Vol.  diel., 

i^  U"ri^  +  (eTri)3 


or  simply  Efcw  =  \|  finz.,  - 

F  U  K+Cj^  4C2  +. . .  Cn 


For  pulsing,  the  ims  of  the  pulse  alone  (which  determines  breakdown)  is  given 

by 

45.  Efp  =  ]j  Pjp-  Pout  7tj[\ 

||  (K  +  C^+  Cg+. ..Cn)  v  Tp  / 

The  resulting  values  of  Efp  for  the  author's  work  were  converted  to  an  equiva¬ 
lent  EfCw  via  Figure  9  for  plotting  and  comparison  with  theory. 


EVALUATION  OF  CONSTANTS. 

C,  in  equation  45  represents  the  Pyrex  vacuum  system  within  the  cavity. 
Because  of  the  power  limits  of  the  equipment,  no  breakdown  could  be  initiated 
at  atmospheric  pressure  and  this  meant  that  K  could  not  be  independently  deter¬ 
mined.  Instead,  C  and  K  were  lumped  and  the  sum  was  found  empirically  to  equal 
1.65  x  10  (corrected  for  pulsing).  C  equals  1.26  x  10“3  from  equation  4  ., 
giving  a  K  of  0.39  x  10 “3  by  subtraction,  versus  a  theoretical  K  of  0.37  x  10-3 
by  equations  37 ,  3^  and  39  as  stated  previously.  Thus  we 
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have  an  expression  for  E  as  a  function  of  measured  power  and  some  constants. 
■The  power  is  read  simply  by  couplers  or  meters,  and  once  the  constants  K,  Cj  , 
Cg,.,.Cn  are  known  a  simple  slide  rule  calculation  is  all  that  is  needed. 

The  (P  in  -  P  out)  term  was  found  by  multiplying  calorimetricaUy  measured 
forward  power  by  1  -(Power  reflected/Power  forward). 

The  ratio  (Prefl/  P  for)  is  found  from  the  oscilloscope  display  of 

rectified  rf  from  the  bidirectional  coupler  as 

i  6  P  refl  3-  (Voltage  reflected)3 
P  for  (Voltage  forward)1 

The  voltage  readings  are  taken  just  before  breakdown,  because  the  corona 
increases  power  absorption.  Figure  10  shows  characteristic  oscilloscope  traces 
of  rf  power  sampled  by  the  bidirectional  coupler.  10a  shows  forward  power  and 
10b  shows  the  relative  reflected  power  (about  25%)  from  the  cavity  containing 
the  pyrex  vacuum  system.  10c  shows  the  reflected  power  during  gas  breakdown 
(not  on  same  scale  as  10a  and  10b  -  the  breakdown  reduces  reflected  power). 

The  actual  breakdown  points  were  determined  by  holding  either  pressure 
or  power  constant  and  varying  the  other  until  a  glow  was  visible  in  the  cavity. 
The  vacuum  system  was  flushed  after  each  breakdown  and  allowed  to  cool  again 
to  ambient  temperature  (besides  dielectric  heating  of  the  Pyrex,  the  glowball 
can  cause  intense  heating  by  conduction-convection). 

A  was  estimated  by  two  methods: 

(1)  When  pA>10,  the  curve  of  log  E  vs  log  P  becomes  linear,  and 
Ee/p  =  32.  Thus  A  =  10/p  =  0.67  cm 

(2)  The  diameter  of  the  smallest  possible  "glowball 11  is  about  one 
inch.  Thus  from  equation 

28.  A  =  r/rr  =  cm/3.14  =  0*61  cm. 

The  value  of  A  from  (l)  is  considered  more  accurate  by  the  author  and  was 
used  for  the  theoretical  curves  of  Figure  11  derived  from  the  simplified 
breakdown  theory  of  MacDonald.  (Fig.  5). 


Comparison  of  Experiment  and  Theory 

Figure  11  shows  the  experimental  breakdown  curve  for  air  at  2450  MHz  and 
theoretical  curves  for  air  at  915  MHZ  and  2450  MHz.  Figure  12  shows  the  exper¬ 
imental  breakdown  curves  for  air,  water  and  COg  . 

Mote  that  the  optimum  breakdown 

for  air  occurs  approximately  where  <0  =  vc  (2nf  =  5.3  x  io9p,  or  p  =  2.95  mm) 
and  that  the  theory  works  well  for  "high"  pressures  where  the  assumptions  of 
the  derivation  holds. 

At  "low"  pressure  where  vc«tOthe  assumption  of  a  homogeneous  electric 
field  no  longer  holds,  and  higher  electric  fields  are  needed  than  theory 
predicts  as  the  electrons  are  swept  out  of  the  intense  field  region.  At  the 
inflection  points  of  the  left  ends  of  the  curves  of  Fig  11  &  12  the  oscillation 
amplitude  approximates  the  geometric  A  which  is  larger  than  the  A  representing 
the  intense  field  region.  As  the  pressure  decreases  below  the  inflection  point, 
the  experimental  breakdown  fields  become  smaller  than  theory  predicts,  as  re¬ 
emission  from  the  walls  becomes  more  pronounced.  Below  about  100  microns 
pressure  the  experimental  breakdown  fields  actually  decrease  again. 

This  decrease  at  "low"  pressures  may  be  very  important  near  the  end  of 
the  freeze  dry  cycle  when  the  constants  Cn  for  the  food  are  small  and  the 
pressure  is  also  very  low.  Bleeding  in  a  high  breakdown  strength  gas  and/or 
raising  the  pressure  may  have  great  advantages  during  this  final  stage,  but  a 
tradeoff  must  be  made  between  the  pressure  gradient  for  mass  flow  and  the 
maximum  dielectric  heating  determined  by  the  breakdown  fields.  Fortunately 

the  mass  flow  for  the  last  stage  is  diffusive,  which  depends  upon  the  partial 
pressure  of  water  vapor  rather  than  total  pressure.  Thus  a  higher  pressure  of 
a  second  or' third  gas  should  not  limit  the  mass  flow  until  the  transition  to 
laminar  flow.  Note  in  Figure  12  that  water  vapor  is  a  gas  which  has  higher 
breakdown  fields  than  air  above  2  mm  Hg  but  lower  than  air  below  that  pressure. 
The  higher  region  is  caused  by  the  high  recombination  rate  of  water  vapor  (recall 
that  recombination  becomes  dominant  at  "high"  pressures).  The  low  pressure 

difference  can  be  attributed  to  a  low  collision  frequency  vc. 
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The  agreement  of  the  theoretical  and  experimental  curves  in  Figure  11  is 
also  a  check  upon  the  validity  of  equation  45,  vhich  was  used  to  determine 
experimental  values  of  E  from  power  measurements. 

A  final  test  of  equation  45  was  made  by  placing  a  second  dielectric  load 
Cg  in  the  cavity.  The  block  of  plexiglass  used  had  a  loss  factor  of  .0166, 
about  the  same  as  frozen  beef,  and  a  C  of  0.58  x  10  .  Referring  to  equation 

45,  we  see  that  more  power  is  needed  for  breakdown  at .a  given  pressure.  As  a 
check  on  the  effect  of  the  plexiglass  block  reducing  the  volume  of  the  Pyrex 
flask  and  possibly  affecting  A,  a  block  of  low-loss  polyethylene  of  equivalent 
volume  was  also  tested.  C  for  the  polyethylene  was  .0027  x  10“3 ,  which  is 
neglibible  in  comparison  to  K  +  Cx  =  1.65  x  10-3.  Table  2  summarizes  these 
results  which  show  agreement  between  theory  and  experiment  within  the  limits 
of  experimental  error. 


Table  2  -  Effect  of  Plexiglass  and  Polyethylene  dielectric  Loads  on 
Breakdown  Fields 


Dielectric 
Material  . 

Power 

(watts') 

(Vreflected)3 

"(VforwardJ'3 

Ef  theory  Eq.47 
volts/cm 

Ef  actual 
volts/cm 

Eth/ 

'  Eact 

Plexiglass 

230 

77% 

300 

302 

.995 

e'r  =  2.60 

188 

75% 

267 

260 

I.03 

e"r  =  .0064 

155 

75% 

243 

232 

1.05 

interpolated  between 

1000  MHZ  &  3000  MHZ 

at  25°C. 

128 

75% 

221 

212 

1.04 

Polyethylene 

230 

75% 

339 

338 

1.00 

e"r  =  4  x  10“4 

188 

75% 

307 

309 

.99 

e  !r  —  2.26 

155 

7&% 

279 

284 

.98 

128 

75% 

253 

272 

.93 

110 

75% 

235 

232 

1.01 
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SUMMAiff 


A  method  has  been  derived  and  tested  for  determining  absolute  electric 
field  strengths  and  approbate  breakdown  fields  for  the  general  case  of  freeze 
drying  with  microwaves.  The  parameters  of  frequency,  cavity  shape  and  size, 
pressure,  temperature,  and  dielectric  load  were  included  in  the  analysis.  The 
frequency  has  a  strong  effect  on  the  economics  of  microwave  freeze  drying  and 
2450  MHZ  appears  better  than  915  MHZ  from  a  theoretical  standpoint.  The  use 
of  pulsing  to  avoid  breakdown  is  not  practical  because  of  the  reduced  power 
absorption.  For  precise  prediction  of  maximum  power  input  without  breakdown 
throughout  a  freeze-dry  cycle,  the  gas  temperature,  pressure  and  composition 
must  be  known,  along  with  sublimation  rates  and  dielectric  constants  and  loss 

factors  for  the  food.  The  effect  of  the  partial  pressure  of  water  vapor  is 
important  and  more  work  needs  to  be  done  on  breakdown  in  mixturaof  air  and 
water  vapor.  Future  work  should  also  relate  the  sublimation  rate  to  the 
theoretical  heating  rate  which  is  easily  obtainable  from  equation  44  or 
45  in  conjunction  with  the  power  equation  4l  and  42.  Lastly,  future  work 
should  test  the  derivations  for  food  dielectric  loads  throughout  freeze 
drying  cycles. 

This  methodology  can  be  applied  to  any  other  microwave  heating  application 
as  well  for  determining  electric  field  strengths  and  theoretical  heating  rates. 
The  value  of  K  can  be  calculated  theoretically  by  equation  37,  38  and  39^  (Note 
that  39  refers  only  to  a  TE,oa  mode  -  see  the  reference  21  for  others)  or  by 
experiment.  In  the  latter  case,  operating  at  a  pressure  high  enough  such  that 
P  ^  ^  S'ftd  p  \  >  100,  the  absolute  values  of  the  field  are  simply  Ee/p  =  32 
volts/cm-mm  Hg  and  K  can  be  easily  determined  from  equation  35. 


28 


While  electric  field  strengths  can  be  determined  more  accurately  by 
other  methods  for  waveguides  and  small  cavities,  this  method  is  suitable 
for  large  multi-mode  cavities.  The  main  advantage  of  the  breakdown  method 
over  bolometers,  etc.  is  that  the  electric  field  is  not  disturbed  by  a 
detector. 

Also,  the  use  of  coronas  may  be  useful  in  design  and  testing  of  cavities. 

A  cavity-contained  corona  gives  a  direct  display  of  relative  field  intensities, 
identifying  hot  spots  and  modes  (see  Fig.  1).  This  use  was  first  suggested  by 
the  Amperex  Co.  in  one  of  their  technical  reports35,  but  it  explained  no  means 
of  quantitative  measurements. 

A  final  note  on  the  accuracy  of  the  author's  method:  a  homogeneous  electric 
field  was  assumed  throughout  the  cavity  for  determining  power  losses.  While 
this  is  not  a  true  assumption  for  point  to  point  power  absorption,  taking  the 
total  cavity  absorption  assumes  that  the  maximum  and  minimum  electric  field 
average  out  to  the  same  effect  as  an  average  even  field  throughout  the  cavity. 
This  assumption  proved  to  be  good,  at  least  for  the  cavity  used  in  this  work, 
and  should  hold  for  any  well  designed  cavity  with  a  load  large  enough  to  occupy 
both  maxima  and  minima  in  the  electric  field.  For  applications  to  freeze  drying, 

we  should  also  keep  in  mind  that  the  inaccuracies  due  to  measured  e"r  and  e'r 
for  foods  are  probably  greater  than  the  errors  due  to  assumptions  involved  in 
the  author's  derivations. 

Lastly,  as  described  by  Copson33,  the  dielectric  heating  of  a  Pyrex  vacuum 
system  within  the  cavity  like  that  used  in  this  work  contributes  significant 
radiative  heating  for  longer  cycle  times.  This  makes  such  a  setup  unsuitable 
for  freeze  drying,  and  the  recommended  alternative  is  to  place  the  cavity 
inside  a  larger  vacuum  chamber  as  Hoover  et  al  1  ,2  did  or  by  making  the  chamber 
itself  serve  as  the  vacuum  container. 
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LIST  OF  SYMBOLS 


A  cross  sectional  area 


2 

cm 


a 


2 

a  acceleration  cm/ sec 

b  constant  of  electric  field 
distribution  dimensionless 

D  electron  diffusion  coefficient 
cm  sec”1 

E=Erms  root  mean  square  of 

electric  field  volts/cm 

Eg  effective  electric  field 

volts/cm 

EQ  maximum  electric  field 

volts/cm 

e  electron  charge  coulombs 

f  field  frequency  Hertz 

fp  pulse  repetition  rate 

pulses/sec 

h  efficiency  of  attachment 

attachments/ co 1 1 is ion 

h^  probability  of  ionization 

ioni zat ions/co 1 1 is ion 

hr  probability  of  recombination 

recombinations/coll is ion 

I  current  amperes 

K  admittance  of  cavity  ohm"1 

k  thermal  conductivity 

hr-ftz  F 

m  electron  rest  mass  grams 

N  number  of  molecules 

Ne  number  of  electrons 

N0  Avogadro's  number  6.02  x  10^ 
molecules/mole 

*1, 


y 

« 


A 

A 


Cl) 


0), 


attachment  coefficient  attachments/ 
cm  travel 


ionization  coefficient  ion  pairs 
produced/ cm  travel 
constant  equal  to  D/A^ 

differential  operator  3/3x  +  3/3y  +  3/3z 
half  power  skin  depth  cm 

-14 

free  space  permittivity  8.84  x  10 
farad/ cm 


relative  dielectric  constant 


relative  loss  factor 


phase  angle  between  current  and  voltage 
radian 

characteristic  diffusion  length  cm 

free  space  wavelength  of  applied  field 

attachment  frequency  of  one  electron 
attachments/sec 

collisional  frequency  of  one  electron 
coll is ions/ sec 

ionization  frequency  of  one  electron 
ioni zat ions/ sec 

recombination  rate  of  one  electron 
re comb inat ions/ sec 

applied  field  frequency  radians/sec 

resonant  frequency  radians/sec 


32 


n  electron  concentration 

e  electrons/ cur 

3 

P  power  absorbed  watts/cm 

PQ  power  absorbed  at  resonance 

watts 

p  pressure  mm  Hg 

Q  electron  flux  electrons/sec 

energy  stored  in  cavity 
Qu  ^  energy  lost  in  cavity 

q  heat  flux  Btu/hr 

R  gas  constant 

r  radius  cm 

Tp  total  pulse  cycle  time  sec 

peak  to  peak 

t  time  sec 

tp  pulse  time  per  cycle  sec 

with  non- zero  power 

3 

V  volume  cur 

electron  drift  velocity  cm/sec 

X#Y,Z  dimensions  of  rectangular 
parallel  piped  cm 

Subscripts  p  *  pulsed 

cw  »  continuous  wave 
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FIG.  la.  Microwave  Breakdown  In  Air. 

High  power  showing  maximum  corona  (at  constant  pressure). 


FIG.  lb.  Microwave  Breakdown  in  Air. 

Medium  power  showing  reduced  corona  (at  constant  pressure). 
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FIG.  lc.  Microwave  Breakdown  In  Air.  '  •  '■ . 

Low  power  showing  minimum  corona.  Note  llBO-Mtioh  of  hot  spots 
corresponding  to  maxima  in  electric  field. 


E(  VOLTS /CM) 


FIG.  2.  Microwave  in  nitrogen,  air,  oxygen  hydrogen,  and  argon  at  994  MHz  and 
=  I.51  cm.  A.D.  MacDonald,  D.U.  Gaskell,  H.H.  Gitterman  (1963) 
PHYHEV  J1  V130  P.1841 
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FIG.  3.  Breakdown  fields  as  a  function  of  pressure  for  different  frequencies 
(gas,  Heg,  a  mixture  of  Hie  and  Hg  -  =  0.6  cm).  A.D.  MacDonald, 

Microwave  Breakdown  in  Gases.  Wiley  &  Sons,  NY  P.7. 
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FIG.  5 


Microwave  Breakdown  (Ee  s  p  )  for  Air. 

A.D.  MacDonald,  D.U.  Gaskell,  H.N.  Gittennan  (1963) 
PHYREV.  Jl  V130  P.  l84l 
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FIG.  6.  Experimental  set  up  for  determination  of  microwave  breakdown 
fields  -  block  diagram. 


meter 


PIG.  7*  Photograph  of  experimental  set  up  for  determination  of  microwave 
breakdown  fields. 


FIG.  8.  Microwave  cavity  and  vacuum  system  for  determination  of 
breakdown  fields 
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CYCLES  PER  PULSE 

FIG.  9.  Ratio  of  pulsed  to  CW  breakdown  /for  air  plotted  as  a  function  of  cycles  per  pulse. 

Frequencies  from  0.99  to  24.1  GHz.  Data  are  for  pressures  near  the  minimH.  of  break¬ 
down  curves.  A.D.  MacDonald,  ktcrowave  Breakdown  in  Gases.  Wiley  &  Sons,  NY  P  170 
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Oscilloscope  traces  of  sampled  rf  power 
FIG.  10a.  Forward  power. 
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Oscilloscope  traces  of  sampled  rf  power 
FIG«  10b.  Reflected  power,  no  breaMown 
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Oscilloscope  traces  of  sampled  rf  power 
FIG.  10c.  Reflected  power  -  breakdown. 
47 
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FIG.  11.  Experimental  and  theoretical  microwave  breakdown  in  air, 
=  O.67  cm. 
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FIG.  12.  Experimental  microwave  breakdown  in  air,  vater  vapor,  and 
carbon  dioxide  at  2450  HHz  =  0.6j  cm. 


FOOD  LABORATORY  INTERNAL  DISTRIBUTION  LIST 


Process  Development 


Copies 

25  -  Chief,  Technical  Plans  Office,  ILABS 

(20  for  transmittal  to  Defense  Documentation  Center) 

2  -  Technical  Library,  NLABS 

10  -  Program  Coordination  Office,  Food  Laboratory,  NLABS 

7  -  Division  Chiefs,  Food  Laboratory,  NLABS 

2  -  Marine  Liaison  Officer,  NLABS 

3  -  Air  Force  Liaison  Officer,  NLABS 

1  -  Director,  Earth  Sciences  Laboratory,  NLABS 

2  -  Director,  General  Equipment  and  Packaging  Laboratory,  NLABS 

Director,  Pioneering  Research  Laboratory,  NLABS 
Project  Officer  and  Alternate  Project  Officer,  Food  Lab,  NLABS 


3 


FOOD  LABORATORY  DISTRIBUTION  LIST 


Process  Development 

Copies  Copies 


1  -  Commanding  General 

US  Amy  Medical  Research 
and  Development  Command 
ATTN:  MEDDH-SI 
Washington,  D.C.  20315 

2  -  Commanding  General  1 

US  Army  Test  and  Evaluation  Command 
ATTN:  AMSTE-BC 

Aberdeen  Proving  Ground,  MD  21005 


1  -  Commanding  General  1 

US  Army  Combat  Development  Command 
Combat  Service  Support  Group 
Fort  Lee,  VA  23801 

1  -  Commanding  General 

US  Army  Combat  Development  Command 

ATTN:  CDCMR-0 

Fort  Belvoir,  VA  22060 

1  Commanding  General  1 

US  Amy  Materiel  Command 
ATTN:  AMCRD-J1 
Department  of  the  Amy 
Washington,  DC  20315 


2  -  Commanding  Officer  2 

Edgewocd  Arsenal 
ATTN:  SMUEA -TSTI-TL 
Edgewocd  Arsenal,  MD  21010 


1  -  Commanding  Officer 

US  Amy  Medical  Nutrition  Laboratory 
Fitszimons  General  Hospital 
Denver,  CO  80240 


-  Commander 

Defense  Personnel  Support  Center 
ATTN:  Directorate  of  Subsistence, 
DPSC-STS 

2800  South  Street 
Philadelphia,  PA  19101 

Commandant  of  the  Marine  Corps 
Code  A04D 

Washington,  D.C.  20380 


1  -Commandant  of  the  Marine  Corps 
ATTN:  Code  CDE 
Washington,  DC  20380 


-  Executive  Secretary 
Interdepartmental  Committee  on 

Radiation  Preservation  of  Food 
Consumer  Products  Division  623 
Business  and  Defense  Service  Admin. 

US  Department  of  Commerce 
Washington,  DC  20230 

-  Director 
Development  Center 

Marine  Corps  Development  and 
Education  Command 

ATTN:  Combat  Service  Support  Division 
Quant ico,  VA  22134 


Copies 


Copies 


1  -  Director 

Division  of  Biology  and  Medicine 
US  Atomic  Energy  Commission 
Washington,  D.C.  20545 


1  -  Director 

US  Amy  Advanced  Material 
Concepts  Agency 
Washington,  D.C.  20315 

1  -  Chief,  Life  Sciences  Division 
Amy  Research  Office 
Office  of  Chief  of  Research 
and  Development 
Washington,  D.C.  20310 

3  -  Office  of  the  Coordinator  of  Research 
University  of  Rhode  Is  land 
Kingston,  Rhode  Island  02681 


1  -  Dr.  Herbert  E.  Hull 

Chief,  Food  Microbiology 
National  Center  for  Urban  and 
Industrial  Health 
Food  irotection  Research 
222  East  Central  Parkway 
Cincinnati,  Ohio  45202 


1  -  Library  Southern  Utilization 

Research  and  Development  Division 
Agricultural  Research  Service 
US  Department  of  Agriculture 
P  0  Box  19687 
New  Orleans,  IA  70119 

2  -  Quartermaster  School  Library 

US  Amy  Quartermaster  School 
Fort  Lee,  VA  23801 


2  -  Technical  Library 

USACDC  Institute  of  land  Combat 
301  Taylor  Drive 
Alexandria,  VA  22314 


1  -  US  Department  of  Agriculture 

Division  of  Acquisitions 
National  Agriculture  Library 
Washington,  D.C.  20250 

2  -  US  Amy  Research  Office 

ATTN :  Technical  Library 
3045  Columbia  Pike 
Arlington,  VA  22040 


Headquarters  12th  Support  Brigade 
ACofS  Services 
ATTN :  Food  Advisor 
Fort  Bragg,  NC  28307 


4  -  Exchange  and  Gift  Division 
Library  of  Congress 
Washington,  D.C.  20540 


1  -  Stimson  Library  2  - 

ATTN :  Documents  Librarian 
US  Amy  Medical  Field  Service  School 
Frooke  Amy  Medical  Center 
Fort  Sam  Houston,  TX  78234 


. 


SeouritycTassiificatlon" 


,r;  /  i  V  DOCUMENT  CONTROL  DATA  -  R  &  D 

'  sv  (Security  classification  oi  tide,  body  oi  abstract  and  indextngarmotation  muatbe  enteredwhen^he^jiveTallreport^i^claa^llie^) 


t.  oaicinating  activity  (Corporate  author) 

u.b.  Army  Natick  Laboratories 
Natick,  Massachusetts  OI76O 


2«.  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2b.  GROUP 


3.  REPORT  TITLE 


MICROWAVE  APPLICATIONS  TO  FREEZE  DEHYDRATION-GASEOUS  BREAKDOWN  VS  ELECTRIC  FIELD 
STRENGTH  ..  s  ... 


4.  DESCRIPTIVE  NOTES  (Type  of  report  and  inclusive  date*) 


0.  AUTHOR(S)  (First  name,  middle  initial,  last  name) 

Jamies  W.  Gould 
Ernest  M. Kenyon 


e.  REPORT  DATE 

October  1970 


7a.  TOTAL  NO.  OF  PAGES 


57 


76.  NO.  OF  REFS 


34 


8a.  CONTRACT  OR  GRANT  NO. 


6.  PROJECT  NO. 

U66 2708D553 


9A.  ORIGINATOR'S  REPORT  NUM’BER(S) 


96.  OTHER  REPORT  NO(S>  (Any  other  number*  that  may  be  assigned 
this  report) 


71  - 


FL 


FL  -  117 


10.  DISTRIBUTION  STATEMENT 

This  document  has  been  approved  for  public  : 
limited.  ■  ; 

,  < 

release  and  sale;  its  distribution  is  un- 

11.  SUPPLEMENTARY  NOTES 

< 

<  •  1 

12.  sponsoring  military  ACTIVITY 

U.S.  Army  Natick  Laboratories 

Natick,  Massachusetts  OI76O 

13.  ABSTRACT 


The  use  of  microwave  energy  can  reduce  freeze -drying  cycles  to  l/2  -  l/lO  of 
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and  related  to  theory  at  2450  MHZ.  Pressures  studied  were  in  the  range  of  10  to  20 
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